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Abstract 

Using the modern equations of state derived from microscopic calculations, we have calculated 
the neutron star structure. For the neutron star, we have obtained a minimum mass about 0.1 M 
which is nearly independent of the equation of state, and a maximum mass between 1.47 M and 
1.98 M which is strongly dependent on the equation of state. It is shown that among the equations 
of state of neutron star matter which we have used, the stiffest one leads to higher maximum mass 
and radius and lower central density. It is seen that the given maximum mass for the Reid-93 
equation of state shows a good consistency with the accurate observations of radio pulsars. We 
have indicated that the thickness of neutron star crust is very small compared to the predicted 
neutron star radius. 
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I. INTRODUCTION 



For identifying of an astrophysical object as a black hole, it is required to know the 
maximum gravitational mass of a neutron star for stability against collapse into a black 
hole. In the other word, it is expected that below a certain maximum mass, degeneracy 
pressure due to the nucleons is sufficient to prevent an object from becoming a black hole 
1|. Therefore, determining the maximum gravitational mass of neutron star is of special 
importance in astrophysics. 

Observationally, the X-ray pulsars and X-ray bursters can offer direct measurement of 
the neutron star mass, but their accuracy is rather poor and due to large errors, measuring 
masses of them are not very useful. Fortunately, the masses of neutron stars have been 
determined with high accuracy using the binary radio pulsars 
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Theoretically, in calculating a maximum mass of neutron star which should be in agree- 
ment with the precise observations, the equation of state of neutron star matter plays a 
crucial role. Theoretical investigations of the neutron star matter indicates that a con- 
siderable uncertainties exist in the behavior of the equation of state, especially at high 
densities. This is due to using the old non phase-shift equivalent nucleon-nucleon potentials 
H). These uncertainties result a significant uncertainty in the maximum mass of neutron 
star. Therefore having a modern equation of state derived from an accurate many-body 
calculation using modern nucleon-nucleon potentials is of particular importance in the mass 
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determination for the neutron star [11|, HJ, H3j • 

Recently, we have obtained the equation of state of neutron star matter using the mi- 
croscopic constrained variational calculations based on the cluster expansion of the energy 
functional Q, ljj. In these calculations, we have employed the modern two-nucleon poten- 



tials such as the new Argonne AVi 8 [20( and charged dependent Reid-93 [21]. This is a fully 
self-consistent technique which does not bring any free parameter into the calculations and 
its results show a good convergence. In this method, for the asymmetric matter calculations, 
a microscopic computation of asymmetry energy is done and therefore its results are more 
accurate with respect to the other methods which use the semi-empirical parabolic approx- 
imation. In fact, using the modern nucleon-nucleon potentials which are explicitly depend 
on the isospin projection (T 2 ), requires a microscopic calculation 22j, |23|, |24j . In this article, 
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we investigate some physical proper ties of neutron star structure using modern equations of 
state of neutron star matter [18|, [19]. 



II. EQUATION OF STATE AND NEUTRON STAR STRUCTURE 

As it is mentioned in the previous section, the equation of state of neutron star matter 
has a key role in determining the neutron star structure, especially its maximum mass. In 
this work, we study the structure of neutron star using the modern microscopic equations of 
state of neutron star matter employing; the AVi8 an d Reid-93 modern two-nucleon potentials 



as well as the AV14 potential 



251 ] . In our calculations, we also consider the effect of three- 



body force using the UV14 + TNI potential in which the effect of three-nucleon interaction 



is included by adding two density dependent terms into the two-body potential [26[]. The 



procedure of our calculations is discussed in references [18|, [19j. Our results for the equation 
of state of neutron star matter are given in Figure HJ It is seen that by inclusion of three- 
nucleon interaction, we find a stiffer equation of state. From Figure [21 it is seen that for 
the Reid-93, AV14 and AVi 8 equations of state, the speed of sound (C s ) is always less than 
the speed of light, even at high densities. In the case of UV14 + TNI equation of state, 
above p = 1.28 fm -3 (e = 29.61 x 10 14 g/cm 3 ), the sound speed exceeds the speed of light. 
However, as we will see, for this equation of state, the neutron star mass reaches a limiting 
value (maximum mass) below this density. 

Using the equation of state of the neutron star matter, we can calculate the neutron 
star mass and radius as a function of central mass density, e c , by numerically integrating 
the general relativistic equation of hydrostatic equilibrium, Tolman-Oppenheimer-Volkoff 
(TOV) equation fl, 

dP G 2 m(r) +A7tr 3 P(r)/c 2 

re 2 

where 

e = p\E(p) + mc 2 ] (2) 

is the mass density, G is the gravitational constant, and 

ATir' 2 t{r')dr' (3) 

has the interpretation of the mass inside radius r. By selecting a central mass density 
e c , under the boundary conditions P(0) = P c , m(0) = 0, we integrate the TOV equation 
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outwards to a radius r = R, at which P vanishes. This yields the neutron star radius R and 
mass M = m(R). 

For densities greater than 0.05 fm -3 , we use our equations of state presented in Figu re [H 
However, at the lower densities, we use the equation of state calculated by Baym et al. [27J , 
since the details of the equation of state at low densities do not affect our results. 

The neutron star gravitational mass (in solar mass M Q units) as a function of central 
mass density, e c , calculated with the AVis, Reid-93, AVu and UV14 + TNI equations of 
state is shown in Figure [3j Our results show that at low densities, the calculated neutron 
star masses exhibit a minimum (~ 0.1 M©) which depends weakly on the equation of state 
of neutron star matter. It can be seen that at high densities, the increasing of mass becomes 
very slow and finally it approaches a limiting value which strongly depends on the equation 
of state. This limiting value of mass is the maximum gravitational mass of neutron star and 
its corresponding central density is the highest possible value for the neutron star central 
density. A star with the higher central density would be unstable against the gravitational 
collapse to a black hole. A comparison between Figure Q] and Figure [3] indicates that there 
is a one-to-one correspondence between the behavior of equation of state of neutron star 
matter, especially at high densities, and the resulting mass of neutron star. The stiffest 
equation of state leads to the higher mass and lower central densities. It should be noted 
that, as mentioned above, the limiting value of neutron star mass with the UV14 + TNI 
equation of state is reached below the density in which the speed of sound exceeds the speed 
of light. 

In Figure HI We have presented the radius of neutron star with the AVis, Reid-93, AV14 
and UV14 + TNI potentials as a function of central mass density. It is seen that as the 
central density increases, the radius decreases very rapidly and then for e c > 7 x 10 14 g/cm 3 , 
it reaches a nearly constant value. This shows that the radius of neutron star is nearly 
independent of the central density. From Figure HI we also observe that the stiffest equation 
of state leads to a relatively larger radius for the neutron star. 

The gravitational mass versus radius (M-R curve) for the AVig, Reid-93, AV14 and 
UV14 + TNI potentials is shown in Figure [5j Our results indicate that for the neutron 
star, there is a minimum gravitational mass about 0.1 M Q which is nearly identical for 
different equations of state, and a maximum mass which is different for these equations of 
state. It is well known that in this mass region, the equilibrium configuration of neutron 
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stars can exist. 

A summary of our results for the properties of maximum mass configuration of neutron 
star predicted for different equations of state are given in Table IIHI From Table 1111} it can 
be seen that the inclusion of the three-nucleon interaction considerably affects the calculated 
properties of neutron star maximum mass configuration. Here, we explicitly see that the 
stiffest equation of state used in our calculations gives the lower central density and higher 
maximum mass and radius for the neutron star. As it is seen in Table 1111} our results 
give the maximum mass of neutron star between 1.47 M Q and 1.98 M . This agrees with 
;he measured range of neutron star masses and the results of other theoretical investigations 
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28|, |29j. Furthermore, our predicted maximum mass of neutron star with the Reid-93 
potential has a good consistency with the mass determined from the accurate observations 
of radio pulsars j^, l], [9}. This is a good confirmation for our equation of state of dense 
matter with the modern nucleon-nucleon potential, Reid-93. If we use the equation of state 
of noninteracting beta-stable matter, a maximum mass 0.7 M Q is obtained. Comparing this 
value with those presented in Table IIHI implies that at high densities, the nucleon-nucleon 
interaction is sufficiently repulsive to produce a considerable shift in the maximum mass of 
neutron star with respect to the noninteracting case. 

For a neutron star, the binding energy is defined as the energy needed to transform it into 
a dispersed configuration. Formally, the binding energy can be obtained from the difference 
between amu mass Ma and gravitational mass, where Ma is given by 28[] 

f R 47rr 2 e(r) 
MA = mA /o [1 - 2m(r)G/rc^ dr ' (4) 

Here, we take as one amu (1.66 x 1CT 24 g). For all values of gravitational masses (greater 
than the predicted minimum mass), our equations of state give positive binding energies for 
neutron stars. This indicates that the neutron stars are always stable against the dispersed 
configuration. 

We have calculated the mass density versus radial coordinate for a neutron star with the 
mass equal to 1.47 M for the Reid-93 equation of state. This is presented in Figure [6] which 
shows that up to r ~ 7.8 km (within the core), the mass density is nearly uniform, and then 
it drops very rapidly (close to the crust). Our results indicate that the crust thickness is 
about 0.43 km which is only 5.2% of our predicted neutron star radius. This implies that 
the core contains the most fraction of neutron star matter. 
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By solving the TOV equation, we have also calculated the properties of neutron star 
with the gravitational mass equal to 1.4 M . Our results for different equations of state 
are presented in Table HT1 Comparing Tables HTT1 and HT1 indicates that for a 1.4 M Q neutron 
star, our used equations of state give higher radius and lower central density with respect to 
those predicted for maximum mass configuration. The differences between these two cases 
become more significant for the stiffest equation of state. 

III. SUMMARY AND CONCLUSION 

In recent years, a wide range of information on neutron stars are available from X-ray 
pulsars, X-ray bursters and radio pulsars. Therefore, the study of neutron stars structure 
on the basis of the equation of state of dense matter is of special interests in astrophysics. 
Theoretical investigation of neutron star structure is also important, since the observational 
results leads to the constraints on the equation of state of dense matter. 

In present paper, we have computed some properties of neutron star structure using the 
modern equations of state of neutron star matter obtained by microscopic constrained vari- 
ational calculations based on the modern potentials. A minimum value for the gravitational 
mass of neutron star is predicted near 0.1 M Q which is nearly identical for different equations 
of state. It is seen that a higher maximum mass is obtained, when we consider the effect of 
three-nucleon interaction in our calculations. We have shown that the radius of neutron star 
becomes nearly constant for all central densities approximately greater than 7 x 10 14 g/cm 3 . 
The properties of maximum mass configuration of neutron star are predicted. These prop- 
erties change considerably by considering the three-nucleon interaction. Our results show 
that the stiffest equation of state gives the higher maximum mass and radius and lower 
central density for the neutron star. The obtained maximum mass of neutron star is be- 
tween 1.47 M and 1.98 M Q which agrees with the masses determined from observations. 
In addition, the predicted value 1.47 M by our used Reid-93 equation of state is in a good 
agreement with the accurately determined mass of radio pulsars. We have seen that in our 
calculated mass region, the binding energy of neutron star is positive and therefore it is 
stable with respect to the dispersed configuration. We have computed the mass density 
profile for a neutron star with the mass equal to 1.47 M Q for the Reid-93 equation of state. 
It is shown that the neutron star crust has a thickness about 0.43 km which is very small 
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with respect to the obtained neutron star radius. We have also computed the properties 
of a 1.4 M neutron star with the different equations of state. A higher radius and lower 
central density are predicted for a 1.4 M neutron star compared to the values calculated 
for the maximum mass configuration. Finally, a good agreement between our results for the 
neutron star structure and those of other theoretical calculations is observed. 
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TABLE I: Maximum gravitational mass (M max ), corresponding radius (R), central mass density 
(e c ), and central number density (p c ) obtained with the AVis, Reid-93, AV14 and UV14 + TNI 
potentials. 



Potential 


M max (M ) R (km) e c 


(10 14 g/cm 3 ) p c 


(fm : 


Reid-93 


1.47 


8.23 


35.39 


1.7 


AV 14 


1.56 


8.28 


33.93 


1.6 


AV 18 


1.65 


8.79 


31.16 


1.45 


UVu + TNI 


1.98 


9.81 


27.17 


1.2 



TABLE II: Properties of neutron star with the gravitational mass equal to 1.4 M for different 
equations of state. 

Potential R (km) e c (10 14 g/cm 3 ) p c (fin" 3 ) 

Reid-93 8.51 27.71 1.4 

AV M 8.81 22.83 1.2 

AVis 9.58 18.25 0.95 

UV14 + TNI 11.16 11.32 0.6 
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FIG. 1: Pressure (P) versus mass density (e) for the neutron star matter with the Reid-93 (full 
curve), AV14 (dashed curve), AVis (dotted curve) and UV14 + TNI (heavy dotted curve) potentials. 
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FIG. 2: As Figure [Q but for the sound speed (C s /c) versus number density (p). 
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FIG. 3: As Figure [H but for the gravitational mass of neutron star (in units of the solar mass 
M ) versus central mass density (e c ). 
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FIG. 4: As Figure [H but for the radius versus central mass density (e c ). 
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FIG. 5: As Figure [H but for the mass-radius relation. 
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FIG. 6: Mass density as a function of radial coordinate for a 1.47 M neutron star with the 
Reid-93 potential. 
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